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Abstract 

Double beta decay of ^^''Nd and ^^^Nd to the excited states of daughter nuclei have been studied 
using a 400 cm^ low-background HPGe detector and an external source consisting of 3046 g of 
natural Nd203 powder. The half- life for the two-neutrino double beta decay of ^^'^Nd to the 
excited 0^^ state in ^^^Sm is measured to be T1/2 = [l.33toil{stat)toil{syst)] ■ 10^° y. For other 
(Oz^ + 2u) transitions to the 2f, 2^, 2^ , and 0^ levels in ^^''Sm, limits are obtained at the level of 
~ (2 — 8) • 10^° y. In the case of -"^^^Nd only hmits for the (Oi^ + 21/) transitions to the 2'f , 0^ , and 
2^ excited states in ^^'^Sm were obtained and are at the level of ~ (4 — 8) ■ lO^'' y. 
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* Electronic address: 



barabash@itep.ru 



I. INTRODUCTION 



The experiments with solar, atmospheric, reactor and accelerator neutrinos have provided 
compelling evidence for the existence of neutrino oscillations driven by nonzero neutrino 
masses and neutrino mixing (see recent reviews and references therein). These 

results are impressive proof that neutrinos have a nonzero mass. However, the experiments 
studying neutrino oscillations are not sensitive to the nature of the neutrino mass (Dirac or 
Majorana) and provide no information on the absolute scale of the neutrino masses, since 
such experiments are sensitive only to the difference of the masses, Am^. The detection 
and study of neutrinoless double beta (Oiyd(3) decay may clarify the following problems of 
neutrino physics (see discussions in ^, la, l6|: (i) neutrino nature: whether the neutrino is a 
Dirac or a Majorana particle, (ii) absolute neutrino mass scale (a measurement or a limit on 
mi), (iii) the type of neutrino mass hierarchy (normal, inverted, or quasidegenerate), (iv) 
CP violation in the lepton sector (measurement of the Majorana CP- violating phases). 

Double beta decay with the emission of two neutrinos (2z//5/?) is an allowed process of 
second order in the Standard Model. The 2zy/3/5 decays provide the possibility of an exper- 
imental determination of the nuclear matrix elements (NME) involved in the double beta 
decay processes. This leads to the development of theoretical schemes for NME calculations 
3oth in connection with the 2u(3j3 decays as well as the Ou/Sf] decays (see, for example, 
3, y, Q])- At present, 2z^/?/5 decay to the ground state of the final daughter nucleus has been 
measured for ten nuclei: ^^Ca, ™Ge, ^^86^62^^ 100^^^ ii6q^^ i28rjg^ laor^g^ iso^^ ^^^^ 238u 

(a review of the results is given in Refs. Q, [ll|, 12|). 



The f3(3 decay can proceed through transitions to the ground state as well as to vari- 
ous excited states of the daughter nucleus. Studies of the latter transitions allow one to 
obtain supplementary information about (3(3 decay. Because of smaller transition energies, 
the probabilities for /3/3-decay transitions to excited states are substantially suppressed in 
comparison with transitions to the ground state. But as it was shown in Ref. jlSt], by us- 
ing low-background facilities utilizing High Purity Germanium (HPGe) detectors, the 2v(3(3 
decay to the 0^ level in the daughter nucleus may be detected for such nuclei as ^°°Mo, 
^^Zr, and ^^°Nd. For these isotopes the energies involved in the (3(3 transitions are large 
(1903, 2202, and 2627 keV, respectively), and the expected half-lives are of the order of 
]^q20 _ ]^g2i y rjnj^g required sensitivity was reached for ^''^Mo in four independent experi- 
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ments 



It was also obtained in ^^°Nd |19|]. Recently additional isotopes 



('82 c 



Se, ^°Te, Cd, and Ge) have also become of interest to studies of the 2z//3/3 decay to 



the Oi level (see reviews in Refs 
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The OiyjSfJ transition to excited states of daughter nuclei provide a clear-cut signature of 
such decays, and is worthy of a special note here. In addition to two electrons with fixed 
total energy, one (0"'' — 2+ transition) or two (0"'" 0^ transition) photons appear, with 
their energies being strictly fixed. In a hypothetical experiment detecting all decay products 
with a high efficiency and a high energy resolution, the background can be reduced to nearly 
zero. This zero background idea will be the goal^future experiments featuring a large mass 
of the (3P sample (Refs. [20 
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24| it was mentioned that detection of this 



123|). In Ref. 

transition will give us the additional possibility of distinguishing the OujSf] mechanisms (the 
light and heavy Majorana neutrino exchange mechanisms, the trilinear R-parity breaking 
mechanisms etc.). So the search for f3p transitions to the excited states has its own special 
interest . 

In this article, results of an experimental investigation of the (3/3 decay of ^^°Nd and 
^^^Nd to the excited states in ^^°Sm and ^^^Sm are presented. The decay schemes for the 



triplets i^°Nd - ^^°Pm - ^^^Sm [25i] and ^^^Nd - ^^^Pm - ^^^Sm |26j are shown in Fig. [Hand 



Fig. [2] respectively. The search for f3p transitions of ^^°Nd and ^"^^Nd to excited states has 
been carried out using a HPGe detector to look for 7-ray lines corresponding to their decay 
schemes. A preliminary result for (]j3 decay of ^^°Nd to the Oj'' excited state of ^^''Sm was 
published in Ref. |l9|. 



II. EXPERIMENTAL STUDY 



The experimental work was performed in the Modane Underground Laboratory (depth 
of 4800 m w.e.). A 400 cm'^ low-background HPGe detector was used to measure a 3046 g 
sample of Nd203 powder in a special Marinelli delrin box which was placed on the detector 
endcap. Taking into account the natural abundance there are 153 g of ^^°Nd (5.64%) and 
154 g of ^"^^Nd (5.76%) in the sample. Data were collected for 11320.5 h. 

The Ge spectrometer is composed of a p-type crystal. The cryostat, endcap, and the 
other mechanical parts have been made of a very pure Al-Si alloy. The cryostat has a J-type 
geometry to shield the crystal from radioactive impurities in the dewar. The passive shielding 
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FIG. 1: Decay scheme of ""^Nd are taken from 



1255.40 

1193.73 

1046.04 

773.24 

740.38 

333.86 



25l ]. Only the investigated levels of 0"*" and 2"^ and 



levels associated with transitions to the first ones are shown. Relative branching ratios from each 
level are presented. 
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FIG. 2: Decay scheme of ^^^Nd Only the investigated levels of 0^ and 2+ are shown. The 
relative branching from each level are also given. 



consisted of 2 cm of archeological lead, 10 cm of OFHC copper, and 15 cm of ordinary lead. 
To remove ^^^Rn gas, one of the main sources of the background, a special effort was made 
to minimize the free space near the detector. In addition, the passive shielding was enclosed 
in an aluminum box flushed with high-purity nitrogen. 

The electronics consisted of currently available spectrometric amplifiers and a 8192 chan- 
nel ADC. The energy calibration was adjusted to cover the energy range from 50 keV to 3.5 
MeV, and the energy resolution was 2.0 keV for the 1332-keV line of ^°Co. The electronics 
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FIG. 3: Energy spectrum with natural Nd203 powder of observed 7-rays in the range [280-380] 
keV. 



were stable during the experiment due to the constant conditions in the laboratory (tem- 
perature of ~ 23° C, hygrometric degree of ~ 50%). A daily check on the apparatus assured 
that the counting rate was statistically constant. 

The current data of accepted values for different isotopes published in Nuclear Data 
Sheets were used for analysis of the energy spectrum. The photon detection efficiency for 
each investigated process has been calculated with the CERN Monte Carlo code GEANT 
3.21. Special calibration measurements with radioactive sources and powders containing 
well-known ^^^Ra activities confirmed that the accuracy of these efficiencies is about 10%. 

The dominate detector backgrounds come from natural "^^K, radioactive chains of ^^^Th 
and 235,238u^ man-made and/or cosmogenic activities of ^^"^Cs and ^°Co. The sample was 
found to have a large activity of (46.3 mBq/kg). Additionally long-lived radioactive 
impurities were observed in the sample, but with much weaker activities, i.e. ^^''Cs (0.089 
mBq/kg), ^^^Lu (0.450 mBq/kg), ^^^La (0.068 mBq/kg), ^^^B& (0.155 mBq/kg), etc. In our 
case the most important isotopes contributing to energy ranges of the investigated transitions 
are ^^^Bi (1.15 mBq/kg), ^^'^Ac (0.93 mBq/kg), ^^''Ac (0.62 mBq/kg), and their daughters. 

Figures [3l m [5l and [6] show the energy spectra in the ranges of interest. 
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FIG. 4: Energy spectrum with natural Nd203 powder of observed 7-rays in the range [360-460] 
keV. 
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FIG. 5: Energy spectrum with natural Nd203 powder in the ranges of investigated 7-rays ([500-600] 
and [670-770] keV). 



III. ANALYSIS AND RESULTS 

A. Search for processes in ^^°Nd 

Double beta decays of ^'^^Nd to 2+ (333.86 keV), 0^ (740.38 keV), 2+ (1046.04 keV), 2+ 
(1193.73 keV) and 0^ (1255.40 keV) levels in ^^osm have been investigated. 
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FIG. 6: Energy spectrum with natural Nd203 powder for 7-rays in the ranges [850-950] and [1140- 
1240] keV. 

1. Decay to the excited state 

The transition is accompanied by two 7-rays with energies of 333.9 keV and 406.5 keV. 
The detection photopeak efficiencies are equal to 2.30% at 333.9 keV and 2.29% at 406.5 
keV. Fig. [3] and Fig. IHshow the energy spectrum in the ranges of interest. As one can see 
there is an excess of events above the averaged continuous background at the investigated 
energies. Isotopes of natural radioactivity (^^^Pb, ^^^Bi, ^^''Th, and ^^'^Ac), found in the 
spectrum, have 7-lines near these energies. ^^^Bi contributes to both investigated ranges 
through 7-rays with energies of 333.37 keV (0.065%) and 334.78 keV (0.036%) for the 333.9- 
keV peak, and 405.72 keV (0.169%) for the 406.5-keV peak. Ac touches the 333.9-keV 
peak range with its 7 (332.37 keV, 0.40%). Appropriate 7-rays from ^^^Th and ^^^Pb are 
334.37 keV (1.54%) and 404.853 keV (3.78%), respectively. Both of the mentioned isotopes 
are daughters of ^^'^Ac (T1/2 = 21.772 y; ^^^U chain). Activity of ^^'^Ac was estimated using 
the most intensive 7-lines of its daughters, ^^'^Th (235.96 keV, 12.9%; 256.23 keV, 7.0%), 
223Ra (269.46 keV, 13.9%) and ^i^Rn (271.23 keV, 10.8%; 401.81 keV, 6.6%). 

There is also the cosmogenic isotope, ^^°Eu (T1/2 = 36.9 y), which decays to ^^°Sm, with 
7-rays of 333.9 keV (96%), 406.5 keV (0.14%), 439.4 keV (80%), 584.3 keV (52.6%). The 
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TABLE I: Analysis of events in the energy ranges of the peaks under study. 



Peak (keV) 333.9 406.5 

Number of events 779 603 

Continuous background 656.4 it 1.6 484.6 ± 1.2 

Isotope 2i4g-^333 37.334 73) 7.9 ± 2.2 2143^(405 72) 9.3 ±2.1 

contributions 227^^^(334 37) 30.7 ±4.7 2iip|3(404_85) 10.2 ± 1.1 

228Ac(332.37) 5.6 ± 0.7 

Excess of events 78.5 ± 28.4 99.0 ± 24.7 



line at 439.4 keV is within one standard deviation of the continuous background, therefore 
it can be taken into account systematic error. 

Table [I] presents the results of the analysis for the two peak energy ranges being studied. 
A peak's shape is described by a gaussian with a standard deviation of 0.58 keV at 333.9 keV 
and 0.61 keV at 406.5 keV. For the analysis a peak's range is taken to within approximately 
four standard deviations {E±2a), i.e., 94.82% of full peak at 333.9 keV and 92.75% of full 
peak at 406.5 keV. As one can see there is an excess of events for each peak. Summing the 
two peaks we obtain a signal of (177.5 ± 37.6) events, corresponding to a half-life of ^^''Nd 
to the first 0+ excited state of ^^°Sm given by T1/2 = [1.33l:°oil{stat)tojUsyst)] ■ lO^o y. 

The primary systematics come from the GEANT calculations (10%), continuous back- 
ground estimation (2.6%), and the possible contribution of ^^''Eu in the 333.9-keV peak 
(8.5%). 

Previous experiments gave only limits on this transition, > 1 ■ 10^° y and > 1.5 ■ lO^^ 
y 28|. Taking into account all errors, our "positive" result is not in contradiction with the 



previous limits. 



2. Decay to the 2~[ excited state 

To search for this transition one has to look for a 7-ray with an energy of 333.9 keV. The 
detection efficiency is 2.60%. The analysis given above shows that the excess of events at 
333.9 keV is mainly due to the double beta decay of ^^°Nd to the Of excited state of ^^°Sm. 
So one can only give the lower half-life limit on the transition to the 2^ excited state of 
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TABLE II: Experimental results for {0u + 2u)(3(3 decay of ^^°Nd to the excited states of ^^°Sm. All 
limits are given at the 90% C.L. 



Excited state 


Energy of 7-rays 


(ry+2-)e., (1020 y) 






(efficiency) 


this work 


other works 


(333.86) 


333.9 (2.60%) 


> 2.2 


> 0.91[27j 

> 24*^ \31] 


0t(740.38) 


333.9 (2.30%) 
406.5 (2.29%) 


= [l.33+'^nf^(stat)+°fJsyst)]'' 


> 1.5 [28] 

> 1 [27] 

> 2.4'^[3]J 


2+ (1046.04) 


712.2 (1.78%) 


> 8.0 


> 1.4 [32] 


2+(1193.73) 


1193.7 (0.95%) 


> 5.4 


> 0.027 [33] 


0^(1255.40) 


921.5 (1.45%) 


> 4.7 


> 2 [27] 



"Only Oi/ decay mode 

''Half-life value for 2v decay (see text for the details) 



29 



30 



^^°Sm. The limit has been calculated using the likelihood function described in Refs. 
which takes into account all the peaks identified above as background. This result together 
with available data on decay of ^^°Nd from other experimental works are presented in 
Table II. 



3. Decays to the 2^ , 2^ and 0^ excited states 

To search for these transitions one has to look for 7-rays with energies of 712.2, 921.5 and 
1193.7 keV (Fig.[T]). As one can see from figures [5] and [6l there are no statistically significant 
peaks at these energies. Using the same technique as above {29! . I30I the lower half-life limits 
are found within (4.7 — 8.0) • 10^° y for the transitions (Table HTj) . Table HTl also presents 
other valuable data on these transitions. 
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TABLE III: Experimental results for {Ou + 2v)[5[5 decay of ^^^Nd to the excited states of ""^^^Sm. 
All limits are given at the 90% C.L. (See text for details.) 





J— iiici gj' y 1 dij o 

(efficiency) 


/rp0v+2v\ 

^ 1/2 '^^P 

this work 


other works 


2|(550.26) 


550.3 (2.36%) 


> 6.6 


> 0.03 [33] 


0]*" (1424.46) 


550.3 (2.16%) 


> 7.9 






874.2 (1.83%) 






2+(1454.12) 


550.3 (1.11%) 


> 3.8 


> 0.027[33] 




903.9 (0.87%) 







B. Search for (3(3 processes in ^Nd 

A search for the double beta decays of ^^^Nd to the 21, , and 2^ excited states of 
^^^Sm was carried out by looking for 7-rays with energies of 550.3, 874.2, and 903.9 keV 
accompanying these transitions (Fig. [2]). Figures [5] and [6] show no statistically significant 
peaks at these energies. The lower half-life limits reported in Table HIT] have been calculated 
using the same procedure as in section IIII A2[ Available data on 13(3 decay of ^^^Nd from 
other experimental works are also presented in Table lllli 



IV. DISCUSSION 



Because the technique used in the present work does not allow for a distinction between 
Qh'(3[3 and 2i'(3(3 decay, our result for double beta decay of ^^°Nd to the excited state in 
^^°Sm is the sum of the Oz//3/3 and 2^(3(3 processes. However we believe that we detected 
only the 2^(3(3 decay. This conclusion is supported by two arguments. First, in the recent 
NEMO paper 3l(| the limit on Oz^/3/3 decay of ^^°Nd to the excited state was established as 
2.4-10^'^ y, which is stronger than the half-life value obtained here. Second, the experimental 
limit for the Oz^/3/9 decay of ^^°Nd to the ground O"*" state of ^^''Sm is about two orders of 
magnitude larger 311] than the value reported here. Therefore, considering the reduced 
phase space factors available for the transition to the excited 0^ state, it is safe to assume 
that our result for T112 refers solely to the 2z//3/3 decay. As to the possible contribution of 
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{0u+2iy)(3P decay to the 2f excited state to the peak at 333.9 keV one can conclude that this 
contribution is small. First, the experimental limit on neutrinoless mode (T1/2 > 2.4 • 10^^ 



y 



3l|) is much stronger than the obtained "positive" result for the decay of ^^°Nd to the 
0^ excited state of ^^''Sm. And second, the theoretical prediction for two neutrino mode 
(Ti/2 ~ 7.2 • 10^^ — 1.2 • 10^^ y [34]) is very far from the observed value. 

Using the phase space factor value G = 1.2 • 10'^'^y~^ (for g^ = 1.254) [35j and the 
measured half-life T1/2 = [1.33l[!;|(stat)t[j:?^(syst)] • 10^° y, one obtains a NME value for the 
2z//3/3 transition to the Oj/" excited state of M2u{0^) = 0.025 ± 0.003 (scaled by the electron 
rest mass). One can compare this value with the NME value for the 2z^/3/3 transition to 
the ground state of ^^°Sm, M2,^(0t^ = 0.0331q;qq2 (here we used the average half-life value 
Ti/2 = (7.8 ± 0.7) • 10^8 y from ll| and G = 1.2 ■ 10~^^y~^ (for g^ = 1.254) from [Ssj). 
One can see that M2jy(0+^ ) is ~ 25% greater than M2u{0i). Nevertheless, these values are 
very close and taking into account errors and possible uncertainties their equality can not 
be excluded. From the viewpoint of the theory it is important to confirm this difference or 
to rule it out. Future more precise measurements for both transitions and especially for the 
transition to the excited state are needed. 

For double beta decay of ^^°Nd to the 2'^ , 2^, and 0^ excited states of ^^'^Sm the limits 



obtained are ~ 2-5 times better than the best previous results [27|, [S^, |33|] . The limit obtained 
for the transition to the 2^ excited state is 200 times better than the previous limit 33 1. 

The ^^^Nd limit for the transition to the Of excited state was obtained for the first time 
and for transitions to 2f and 23 states a sensitivity was achieved that is ~ 200 times better 
than in Ref. 33 1. 



V. CONCLUSION 



Double beta decay of ^^°Nd and ^^^Nd to the excited states of daughter nuclei was in- 
vestigated with a high level of sensitivity. The half-life for the 2z//?/5 decay of ^^''Nd to the 
excited 0^ state in ^^^Sm is measured to be T1/2 = [1.33tof3{stat)toil{syst)] ■ 10^° y. The 
strongest limits for other transitions were established. The sensitivity of this experiment 
could still be increased by a few times using a pure Nd203 (or Nd) sample. Also further 
increases in the sensitivity could be reached using an enriched Nd sample and a multicrystal 
HPGe installation to study larger masses of Nd samples. 
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